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Abstract. 

We start the paper with a brief presentation of the main characteristics 
of graphene, and of the Dirac theory of massless fermions in 2+1 dimensions 
obtained as the associated low-momentum effective theory, in the absence of 
external fields. We then summarize the main steps needed to obtain the Hall 
conductivity in the effective theory at finite temperature and density, with 
emphasis on its dependence on the phase of the Dirac determinant selected during 
the evaluation of the effective action. Finally, we discuss the behavior, under 
gauge transformations, of the contribution due to the lowest Landau level, and 
interpret gauge transformations as rotations of the corresponding spinors around 
the magnetic field. 



PACS numbers: ll.10.Wx, 02.30.Sa, 73.43.-f 

1. Introduction 

Graphene is a bidimensional array of carbon atoms, packed in a honeycomb crystal 
structure. Actually, each layer of a graphene sample can be viewed as cither an 
individual plane extracted from graphite, or else as an array of unrolled carbon 
nanotubes. Quite unexpectedly, in 2004, stable monolayer samples of such material 
were obtained b and, in 2005, the Hall conductivity was measured in such samples, 
independently, by two groups [2]. More recently, a different behavior of the Hall 
conductivity was reported [3] for bilayer samples. The main difference between the 
behavior of the Hall conductivity of mono- and bilayer samples is in the height of the 
jump around zero carrier density (or, equivalently, chemical potential). 

From a theoretical point of view, the most remarkable feature of graphene is that, 
in a small momentum approximation, the charge carriers or quasi-particles behave as 
two "flavors" (to account for the spin of the elementary constituents) of massless 
relativistic Dirac particles in the two non-equivalent representations of the Clifford 
algebra (corresponding to the two non-equivalent vertices in the first Brillouin zone), 
with an effective "speed of light" about two orders of magnitude smaller than c [4] . 
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In [5] , we showed that a field theory calculation at finite temperature and density, 
based upon ^—function regularization of the Dirac determinant leads, in the zero 
temperature limit, to a sequence of plateaux in the Hall conductivity consistent with 
the measured ones, each time the chemical potential goes through a nonzero Landau 
level. Moreover, it was shown in [6] that two of the three possible combinations of 
phases of the Dirac determinant in both nonequivalent Clifford representations predict 
a behavior around zero chemical potential consistent with the ones measured in mono- 
and bilayer graphene. For a complete presentation of other approaches to the study 
of graphene see, for instance, [7], and references therein. 

This paper presents, in section[2j a brief introduction to the structure of graphene, 
and to the derivation of the continuous Dirac effective theory, in the absence of external 
fields. Section[3] contains a review of our previous results on the subject, with emphasis 
on the role of the phase of the determinant in giving rise to different behaviors around 
zero chemical potential. In section @J entirely new results are presented. In that 
section, we allow for complex chemical potentials, and concentrate on the contribution 
due to the lowest Landau level, in order to study the invariance of the effective action 
under large gauge transformations. By relating gauge transformations to rotations 
in the plane, we analyze the effect of a 2tt rotation for each of the three possible 
combinations of phases in both representations, and identify, in the zero temperature 
limit, the resulting geometrical or Berry's phases. Our conclusions are presented in 
section [5j 



2. Structure of graphene. Effective continuous model 

In this section, we sketch the main steps leading to the effective Dirac model for 
graphene, in the absence of external fields. For more detailed presentations see, for 
instance, [4]. 

The structure of the direct lattice for graphene is shown in figure [TJ The direct 
lattice is a superposition of two triangular lattices, A and B. The generators of lattice 
A are ai = y/ZaQ^, —^ L ), and a 2 ~ \fZa{^, -^), where a is the lattice spacing. The 

vectors Si = a(0, — 1), s 2 = a(^, \) and s 3 = a(^j^, i) connect each site in the 
lattice A to its nearest neighbor sites in the lattice B. The tight binding Hamiltonian 
can then be written as 

3 

Ho = -t^2Y, i a Hr)b(r + s,-) + tf(r + Si )a(r)] , 

where t is the uniform hopping constant. 

In momentum space, with i? = e~ lkr <?f r ]l, it reads 

l & (k)J r fx A \K r )J 

H = J2 ( $ ( k ) flt ( k ) b ( k ) + $ * ( k ) 6t (k) a (k)) 

k 

3 

with $ (k) = — t ^2 e lk ' Si . After defining two-component spinors as 

i=i 

V(k) = ( a (k), 6(k)) T ,^t(k) = ( a t( k ), &t(k)), one gets 

/ $( k ) 
I $*(k) o 
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Figure 1. Direct lattice for graphene 

Hq vanishes at the six corners of the first Brillouin zone. Among these, only two 
are inequivalent, and can be chosen as 

k=K±=± G^b'°) ; ' t<K±, = a 

When Ho is linearized around these two points , k = K± + p, one obtains 

This last expression shows that each Fermi point gives rise to an effective Dirac 
theory, with an effective "velocity of light" c = ^n, in one of the two inequivalent 
representations of the gamma matrices. Thus, the total Hamiltonian can be taken as 
the direct sum of both (equivalently, as the Dirac Hamiltonian in the reducible 4x4 
representation of the Clifford algebra). Moreover, an overall factor of two (two fermion 
species or "flavors") must be included to take the spin of the original particles into 
account. 

3. The Hall conductivity from the interacting quantum field theory at 
finite temperature and density 

As shown in our previous work on the subject [5J [5], the Hall conductivity can 
be determined by first evaluating the partition function (equivalently, the effective 
action) for massless Dirac fermions at finite temperature and density, in two spacial 
dimensions, in the presence of an external magnetic field perpendicular to the plane, 
and then performing a boost to a reference frame with orthogonal electric and magnetic 
fields. In this section, we sketch our main results in those references, with emphasis on 
the role played by the phase of the Dirac determinant, which appears when treating 
the infinite tower of states associated to the lowest Landau level. We first consider a 
single flavor, and one of the two nonequivalent representations of the Clifford algebra. 

In order to consider the effects due to finite temperature and density, we study the 
theory in Euclidean three-dimensional space, with a compact Euclidean "time" < 
< P, where (3 = (here, ks is the Boltzmann constant and T is the temperature). 
We introduce the (real) chemical potential and the magnetic field through a minimal 
coupling of the theory to an electromagnetic potential = (— i^, 0, Bx\). Natural 
units (c = h = 1) will be used, unless otherwise stated. 

In this scenario, the Euclidean effective action is given by logZ = log det(if) — 
cA)ap, where the subindex AP indicates that antiperiodic boundary conditions must 
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be imposed in the x direction, in order to ensure Fermi statistics. Now, this is a 
formal expression, which we will define through a zeta- function regularization [S], i.e., 



S eff =lo g Z^-± 



C{S, = --T 

s=0 a ds 



where lo represents the eigenvalues of the Dirac operator acting on antiperiodic, 
square-integrable functions, and a is a parameter introduced to render the zeta 
function dimensionless (as expected on physical grounds, our final predictions will 
be a-independent). 

So, in order order to evaluate the partition function, we first determine the 
eigenfunctions, and the corresponding eigenvalues, of the Dirac operator. We propose 

y ktl (x ,x 1 ,x<t) = 7^=8 ( v . Aj = (2Z + l)-. 

V27T/3 V Xk,l[Xl) J p 

Note that, in the last expression, A;, I = — oo, oo are the Matsubara frequencies 
adequate to the required antiperiodic conditions, while the continuous index k 
represents an infinite degeneracy in the direction. 

The resulting spectrum has two pieces: An asymmetric piece, associated to the 
lowest Landau level of the Hamiltonian: 

7T 

uj[ = Xi, with A/ = (21 + 1)— + ifj, and I = — oo, oo , 

P 

with corresponding eigenfunctions 

. , . / (eB)i e -¥(-i-^) 2 \ 

iPk,l{xi) = I q ) ' ^ 

and a symmetric piece 



LJi >n = ± y Xf + 2neB with n=l,...,oo I = — oo, oo , 

corresponding to eigenfunctions with both components different from zero. In all cases, 
the degeneracy per unit area is given by the well known Landau factor, Al = Sp. 

The asymmetric part of the spectrum is quite particular. In the first place, as seen 
from the corresponding eigenfunction is an eigenfunction of the Pauli matrix 03, 
with eigenvalue +1. The eigenfunction with the opposite "chirality" was eliminated 
by the square integrability condition in x\. Moreover, in the other nonequivalent 
representation, this part of the spectrum appears with the opposite sign. Such 
transformation is equivalent to [i — * — /i. This is nothing but charge conjugation 
(which, for real fx, is also equivalent to a parity transformation [S]). As we will discuss 
in what follows, this part of the spectrum is the one which requires the consideration 
of a phase of the determinant when evaluating the effective action. 

Before going to such evaluation, it is interesting to note the invariance of the 
whole spectrum under /x — > /i + This invariance is a natural one, since such 

transformations preserve the antiperiodicity of the eigenfunctions and, thus, the Dirac 
statistics. They are nothing but the so-called large gauge transformations. The 
invariance of the effective action under such transformations is also required for 
topological reasons. We will discuss this point in more detail in section [4] 

As is clear from ([2]), in evaluating the effective action, one must perform the 
analytic extension of the contributions to the zeta function coming from both pieces 
of the spectrum, 



I — — OO 



(2 | + i)-L +i £ 

ap a 
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and 



C 2 ( S .ji,eB) = (l + (-l)- s )A L 

n = 1 

I = — oo 



2ne5 



+ ((2/ + l)- + 



The analytic extension of C2(s, M,ei3) is quite standard, and it relies mainly on 
performing a Mellin transform and making use of the inversion properties of the Jacobi 
theta functions. A detailed presentation can be found in [5]. The final result for the 
contribution to the effective action coming from this piece is (always considering only 
one representation of the gamma matrices and one fermion species) 



Slj f = A L pV2e~B( R 



!-A, > i g[(l + e -(^ 5 ^)' 3 )(l + e-(^i?+^' 3 )] . (4) 



n=l 



The other nonequivalent representation of the Clifford algebra gives rise to an 
identical contribution, since this part of the spectrum is the same in both irreducible 
representations. 

As said before, the extension of (i(s, M> e B) requires a careful consideration of the 
phase of the determinant. In fact, £i can be written as 



Ci(s,M) 



2tt 

a/3 



E 

.1=0 



(I- 



— ) + % — 

T 2tt 



E 

(=0 



(J 



! 2tt 



(5) 



and the definition of the overall minus sign in the second sum depends on the selection 
of the cut in the complex plane of eigenvalues. As discussed in detail in [6], the usual 
prescription is to choose the cut such that one does not go through vanishing arguments 
when continuously transforming eigenvalues with positive real part into eigenvalues 
with negative real part [10] . This prescription then gives rise to what will be called 
in the following the standard phase of the determinant (characterized from now on by 
k = —1). One could certainly choose the opposite prescription, which we will call the 
nonstandard phase (k = +1). Once one of the phases is selected, the contribution of 
Ci to the effective action can be evaluated by making use of the well-known properties 
of the Hurwitz zeta function, to obtain 



Si ff ( K ) = A L 



log 



2cosh(^) 
v 2 ' 



+ K 



When this last contribution is added to the one in (j4|), one gets for the effective 
action 

V 



Seff(K) 



log 



2cosh(^) 



OO 

^2 log "(l + e -(^^B-v)0^ (l + e -(V2^B+^)p^ 



From this last expression, the finite-temperature charge density can be obtained 
as jo (ft) = ~frjjl^ 8.£- I n the zero-temperature (/? — ► oo), and recovering physical 
units, it reduces to 

f(2ec 2 hBn < fi 2 < 2eBc 2 h(n + 1)) = — sign(^) , 
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Figure 2. Hall conductivity for different selections of the phase of the 
determinant. Left to right: K = 1 , K = 2 , K = . In all cases, the horizontal 
axis represents Vq = sgn(/x) fi 2 /2eBhc 2 and the vertical one, a X y h/Ae 2 . 



where n = [ 2 £" Bn ], and [x] is the integer part of x. 

In order to obtain the Hall conductivity, one must perform a boost to a reference 
frame with crossed electric and magnetic fields. The final contribution to the Hall 
conductivity from each fermion species and one irreducible representation is given by 
[6] a xy = -("+— ) £ sign(p) . 

Now, the phases of the determinant in both irreducible representations can be 
selected with the same or with opposite criteria. When this is taken into account, and 
an overall factor of 2 is included, to take both fermion species into account, on obtains 
for the total zero-temperature Hall conductivity 

-4(»+f)e 2 . 
cr xy = j-^ sign{fi) , 

where K — corresponds to selecting the standard phase of the determinant in both 
irreducible representations, K = 1 corresponds to choosing opposite criteria for the 
phases, and K — 2, to choosing both phases in the nonstandard way. The dependence 
of the Hall conductivity on the classical filling factor (vc) is presented in figure^ for 
the three values of K. From that figure, it is clear that the behavior of monolayer 
graphene, as presented in [2], corresponds to K = 1, i.e., to choosing opposite phases 
of the determinant in both representations (or, equivalently, ignoring the phase in both 
representations, as done in [H]). In fact, in this case the (rescaled) Hall conductivity 
shows a jump of height 1 for vc — 0, and further jumps of the same magnitude for 
vc = ±1, ±2, .... In turn, the behavior of bilayer graphene, as reported in [3] is exactly 
reproduced by K = 2 (nonstandard selection of the phase in both representations). 



4. Invariance under large gauge transformations. Interpretation in terms 
of rotations 

To analyze the physical meaning of the invariance of the effective action under large 
gauge transformations in this context, we go back to the zeta function associated to 
the asymmetric part of the spectrum, for two fermion species and one representation, 
this time allowing for an imaginary part in the chemical potential, p, = fi + i 7, while 
always keeping fi ^ 0. In this case, one must be careful when splitting the infinite sum 
as in © (detail of the calculations will appear in T2|). In fact, such splitting must be 
different for different 7-ranges, to make sure that all the eigenvalues in each infinite 
sum have a real part with the same sign, which is crucial in defining the phase. For 
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example, for — A < 22 < I ne has 



<?' ( 1 <r 7/3 <r 1 ^ - 9 A d 




l)7r//3 + i/i- 7 ]" 



E 

Z=0 



e 



[(2/ + 1)tt//3 + + ry) e" 



Now, the values of 9 such that the second term in the RHS does vanish are those ones 
for which, simultaneously, (21 + l)n/{3 + [i sin 9 — 7 cos 9 = = fi cos 9 + 7 sin 0. As 
before, we consider here two different definitions of the phase of the determinant, which 
correspond to the standard definition for the phase k = — 1 , and to the nonstandard 
one k = +1. With each one of these prescriptions, the contribution of the asymmetric 
spectrum to the effective a action is given by 



+ log (e~^ (M+ i 7)(l+sgiiM) _|_ g| (M+i7)(l-sgn/i)^ | _ (g) 

Things are entirely different for ^ = ±-|. In this case, one mode in the infinite 
sum defining the zeta function has a vanishing real part. A careful treatment shows 
that, at such points, S^f is discontinuous. For instance, ^//(if = +5) comc ides 
with lim 7/3 1 - of © . An equally carefully treatment of the case 5^ = — \ shows 
that SljfC^ = — 5) = Slff(?^ = f )• This analysis can be extended to other ranges 
of variation of 3^, to obtain, with fc = — 00, 00, 



- g) < ^ < ( fc + 2)) = 2A l [ 2 Sgn ^ + z ( 7 ~ "X )1 

+ log ^ e ~^ ( ' 1+l(7 ~ 2 F L))(1+sgnAl) + e ^ {lJ - +i{ ' y -^ L)){1 - asn ^A I . (7) 

This expression shows that the contribution to the effective action of the 
nonsymmctric part of the spectrum, in this representation of the gamma matrices, is 
invariant under large gauge transformations, no matter which phase of the determinant 
is selected. As already said, such transformations must constitute an invariance. In 
fact, an increase of 17 in the chemical potential corresponds to the multiplication of the 
eigenfunctions with a phase, i.e., ipk,i(x) — > e I7a:o -0fe^(x). So, an increase 17 = -4p 
is a pure gauge transformation which, moreover, preserves the antiperiodicity in xq. 

Due to the fact that these eigenfunctions are eigenfunctions of 03 such that 
C3ip n =o(x) = ijj n =o(x), one can equivalently write gauge transformations in the form 
i>k,l( x ) — * e l ^~ 2lXo 4>k,i{x)- This last expression shows that, as Xq grows from to /3, 
spinors are rotated by 2"f[3, since ^ is the generator of rotations in the plane XiX2- 
In particular, 7 = =£■ corresponds to a 47r-rotation around the magnetic field. Note 
that not only the partition function, but the Abelian Chern-Simons term (and, thus, 
the effective action) is invariant under large gauge transformations. On the other 
hand, 7 = 5 corresponds to a 27r-rotation. At finite temperature, such transformation 
changes the statistics to a bosonic one. For k = +1, it also gives rise to an overall phase 
of 7r per unit degeneracy in the partition function. Such phase is the contribution which 
survives in the zero temperature limit. The three possible combinations of phases of 
the determinant then give a total phase per unit degeneracy in the partition function 
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of 7r (K — 1, which reproduces the behavior of the Hall conductivity for monolayer 
graphene), or (both for K = and K — 2, this last reproducing the behavior of the 
Hall conductivity for bilayer graphene). At his point, it is interesting to note that, in 
order to have a zero-temperature partition function invariant under rotations of 2tt for 
monolayer graphene, the reduced flux through a unit cell of area must be given by 
JL = fiAi = N, with N a positive integer. This is precisely the condition for physical 
states to transform as unidimcnsional ray representations of the magnetic translation 
group [5]. 

5. Conclusions 

The first conclusion, as already stated in [6j, is that two of the three possible 
combinations of phases give behaviors of the Hall conductivity which coincide with the 
ones measured in mono- and bilayer graphene. In the case of bilayer graphene [3J, the 
(rescaled) Hall conductivity presents a jump of height 2 for vq — 0, and further jumps 
of height 1. The main point here concerns the positions of these subsequent jumps. 
As a matter of fact, according to figure l.b in [3J, these subsequent jumps appear 
for vc = ±1,±2,..., which is exactly the behavior predicted, in our calculation, for 
K = +2. However, the same reference interprets the Hall behavior of bilayer graphene 
through a theoretical prediction first made in [T3J which, as discussed in [2] , predicts 
a plateaux of larger width. Our calculation, instead, completely coincides with the 
measured behavior of the plateaux, both in height and width. 

An entirely new conclusion is that, in each representation, the effective action 
per unit degeneracy is invariant under large gauge transformations, with any of the 
two possible selections of phase. As a result, the invariance persists no matter which 
of the three possible combinations of phases is selected. Moreover, each of the two 
selections of phase in each representation corresponds to a different geometric phase 
under the rotation of spinors along a closed path around the magnetic field (k = — 1: no 
geometric phase; k = +1: geometric phase of tt). So, different values of K correspond 
to different total geometric phases per unit degeneracy, to be compared with the 
Berry phases studied, for instance, in [T3]. Finally, by taking 7/? = -|, we note that, 
for = 1, these three values of K also correspond to the three nonequivalent unitary 
representations of the generator of the cyclic group C3 , which is the relevant symmetry 
in the case of free graphene. 

To the best of our knowledge, the relation between the phase of the fcrmionic 
determinant and Berry's phase hadn't been noticed before. This point, as well as the 
connection with the magnetic translation group will be studied in more detail in [12] . 
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